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Abstract 

Background: Humans can spend the majority of their time indoors, but little is known about the interactions between 
the human and built-environment microbiomes or the forces that drive microbial community assembly in the built 
environment. We sampled 16S rRNA genes from four different surface types throughout a university classroom to 
determine whether bacterial assemblages on each surface were best predicted by routine human interactions or by 
proximity to other surfaces within the classroom. We then analyzed our data with publicly-available datasets representing 
potential source environments. 

Results: Bacterial assemblages from the four surface types, as well as individual taxa, were indicative of different source 
pools related to the type of human contact each surface routinely encounters. Spatial proximity to other surfaces in the 
classroom did not predict community composition. 

Conclusions: Our results indicate that human-associated microbial communities can be transferred to indoor surfaces 
following contact, and that such transmission is possible even when contact is indirect, but that proximity to other 
surfaces in the classroom does not influence community composition. 
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Background 

In the developed world, humans spend a majority of their 
lives indoors. While indoors we encounter microorganisms 
on virtually every surface we touch, and this frequent 
exposure to indoor microbes carries with it the potential 
for disease transmission, as well as interactions with our 
own commensal microbiome [1-3]. Yet we have very lit- 
tle knowledge regarding the ecological processes that 
drive microbial community assembly indoors, nor do we 
understand the degree to which humans share microbes 
with indoor surfaces. 

For any given indoor surface, microbial community 
composition is likely shaped by habitat-specific environ- 
mental constraints (such as the type of surface material), 
and dispersal sources (which include humans, bioaerosols 
and other surfaces within a space). If dispersal among sur- 
faces is a primary determinant of community structure. 
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then adjacent surfaces should be more similar in commu- 
nity composition than surfaces further apart. If dispersal 
from humans is a major determinant, then community 
structure should vary with frequency and nature of human 
contact. 

Because humans harbor distinguishably different micro- 
bial assemblages on different parts of their bodies [4-6], it 
is reasonable to assume that different indoor surface types 
could harbor different microbial communities due to 
frequent contact with specific body parts. Several studies 
suggest that this may be the case. In a recent survey of 
public restroom surfaces, Flores and colleagues [7] found 
that microbes on restroom surfaces were similar to those 
found on specific human body parts, with the strongest 
association observed between toilet surfaces and gut and 
vaginal communities. These associations are to be expected 
given the direct contact common in this environment be- 
tween surfaces and the human body. Associations between 
human use and bacterial community composition have also 
been found on residential kitchen surfaces [8], with bacter- 
ial taxa commonly found on human skin predominating on 
kitchen surfaces, consistent with frequent skin-to-surface 
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contact. It is not known whether more indirect skin-to- 
surface contact, such as sitting fully clothed on a chair, 
might result in similar bacterial transmission. 

Here we test: (a) if different types of surfaces in a 
classroom vary in microbial community composition; 
(b) if their composition varies predictably with the type 
of human contact; and (c) whether these associations 
are stronger than the association between community 
similarity and spatial distance. We sampled desktops, 
floors, walls and chair seats in a classroom to determine 
whether these groups of surfaces harbored distinct bacter- 
ial assemblages related to their differing types of human 
contact. 

Methods 

All samples were collected in the same classroom at the 
Lillis Business Complex, University of Oregon, Eugene, 
OR, USA, during 4-5 August 2011. During the week prior 
to the sampling period, the classroom was occupied by 
classes daily, including on the first day of sampling. Desks, 
floors, chair seats and walls were sampled {n = 18, 18, 
18 and 16, respectively) following identical protocols. 
Samples were distributed throughout the classroom 
(9.2 m X 15.2 m). Pairwise spatial distance between sam- 
ples ranged from 0 to 16 m, with sets of chair, desk and 
floor samples generally taken within 1 m^ of one another. 
Figure 1 displays the spatial distribution of samples within 
the classroom. Surfaces were sampled using a nylon flocked 
swab (copanusa.com; #552C) moistened with sterile buffer 
solution (0.15 M NaCl, 0.1 Tween20). A 289 cm^ (17 cm x 
17 cm) area was swabbed from each surface. Chairs were 
sampled in the center of their upholstered seat surfaces; 
desks (plastic laminate surface) and floors (linoleum) were 
sampled directly above and below chair seats; and walls 
(latex paint) were sampled around the perimeter of the 
room (Figure 1). We were not able to ascertain the 
cleaning schedule prior to sampling, but all surfaces 



were visibly devoid of standing dust. AU samples were 
frozen (-80°C) until DNA extraction. 

DNA extraction, amplification and Illumina library 
preparation followed methods described previously [9,10]. 
DNA was extracted from swabs using a PowerWater 
DNA extraction kit (MoBio Laboratories, Inc., Carlsbad, 
CA, USA) with the following modifications: samples were 
frozen and thawed for two cycles; bead beating length was 
extended to 10 minutes; and samples were eluted in 50 [iL 
Solution PW6. 

The V4 region of the 16S rRNA gene was amplified 
using the F515/R806 primer combination (5'-GTGCCAG 
CMGCCGCGG-3 ', 5 '-TACNVGGGTATCTAATCC-3 '; 
[11]). Amplification proceeded in two steps using a custom 
Illumina preparation protocol where PCRl was performed 
with forward primers that contained partial unique 
barcodes and partial Illumina adapters. The remaining 
ends of the Illumina adapters were attached during 
PCR2, and barcodes were recombined in silico using 
paired-end reads. Adapter sequences are detaUed in 
Additional file 1. AU extracted samples were amplified 
in triplicate for PCRl and triplicates were pooled before 
PCR2. PCRl (25 (iL total volume per reaction) consisted of 
the followdng steps: 5 (iL 5x HF buffer (Thermo Fisher Sci- 
entific, Waltham, MA, USA), 0.5 (iL deoxyribonucleotide 
triphosphates (10 mM, Invitrogen, Life Technologies, 
Grand Island, NY, USA), 0.25 (iL Phusion Hotstart II 
polymerase (0.5 units; Thermo Fisher Scientific), 13.25 (iL 
certified nucleic-acid free water, 0.5 (iL (10 (iM) forward 
primer, 0.5 (iL (10 (iM) reverse primer, and 5 (iL template 
DNA. The PCRl conditions were as follows: initial de- 
naturation for 2 minutes at 98°C; 22 cycles of 20 seconds 
at 98°C, 30 seconds at 50°C and 20 seconds at 72°C; and 
72°C for 2 minutes for final extension. After PCRl, the trip- 
licate reactions were pooled and cleaned with the Qiagen 
MinElute PCR Purification Kit according to the manufac- 
turers protocol (Qiagen, Germantown, MD, USA). Samples 
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were eluted in 11.5 \iL Buffer EB (10 mM Tris-Cl, pH 8.5, 
Qiagen). For PCR2, a single primer pair was used to add 
the remaining Illumina adapter segments to the ends of the 
concentrated amplicons of PCRl. The PCR2 (25 [iL volume 
per reaction) consisted of the same combination of reagents 
that was used in PCRl, along with 5 [xL concentrated PCRl 
product as a template. The PCR2 conditions were as fol- 
lows: 2 minutes denaturation at 98°C; 12 cycles of 20 sec- 
onds at 98°C, 30 seconds at 66°C and 20 seconds at 72°C; 
and 2 minutes at 72°C for final extension. 

Amplicons were size selected by gel electrophoresis, 
extracted, quantified on a Qubit 2.0 Fluorometer (Life 
Technologies, Carlsbad, CA, USA), concentrated, com- 
bined in equimolar concentrations and sequenced on the 
Illumina MiSeq platform at the Dana-Farber/Harvard 
Cancer Center DNA Resource Core (Boston, MA; dnaseq. 
med.harvard.edu). 

Sequence processing was performed using the FastX 
Toolkit (http://hannonlab.cshl.edu/fastx_toolkit) and QIIME 
[12]. Quality filtering settings were as follows: a minimum 
30 quality score over at least 75% of the sequence read, 
no ambiguous bases, and 1 primer mismatch allowed. 
After quality control and barcode assignment, the remaining 
high-quality sequences (mean 13,786 sequences per sam- 
ple ±4,735 SD) were binned into operational taxonomic 
units (OTUs) at a 97% sequence similarity cutoff using 
uclust [13]. High-quality sequences from each OTU cluster 
were taxonomically identified using reference sequences 



from the Greengenes database (2011 release; [14,15]), and 
plant-chloroplast OTUs were removed by name ("Strep- 
tophyta") based on Greengenes taxonomic classifications 
(mean = 8% of sequences in each sample ±10% SD). All 
samples retaining more than 4,000 sequences {n = 15, 14, 
15 and 14 for desks, walls, floors, and chairs, respectively) 
were rarefied to that level for even sampling depth. Se- 
quence files and metadata for all samples used in this study 
have been deposited in Figshare (http://dx.doi.org/10.6084/ 
m9.figshare.687155). Metadata, the unrarefied OTU table, 
and corresponding taxonomic classifications have all been 
included as Additional files 2, 3 and 4, respectively. 

All statistical analyses were performed in R [16], primar- 
ily utilizing multivariate community ecology procedures in 
the labdsv and vegan packages [17,18]. A full record of all 
statistical analysis is included as Additional file 5, and was 
created using the knitr package in R [19]. Community dis- 
tance between samples was calculated using the Canberra 
distance, implemented in the vegan package. Given that 
the four surface types we sampled represent potentially 
similar microbial habitats, we did not expect to see major 
phylogenetic differences among the communities. Thus, 
we used the Canberra taxonomic metric as a way to 
emphasize community differences driven by relatively 
rare taxa with a decreased emphasis on total abundance 
and phylogenetic community distance. Distance-based 
redundancy analysis (DB-RDA) was performed with the 
capscale function in vegan. Individual OTUs exert varying 
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Figure 2 Surfaces harbored significantly different bacterial communities and were linked to differential human contact, (a) Bacterial 
communities were constrained by four different surface types using distance-based redundancy analysis (DB-RDA; constrained inertia = 1 1.4%) 
and were significantly different among types based on Canberra taxonomic distances (P = 0.001 from permutational multivariate analysis of vari- 
ation), (b) Bacterial operational taxonomic units (OTUs) from DB-RDA are shown weighting communities in the same four primary directions. The 
first and second axes from DB-RDA are used in both ordinations (CAP 1 and CAP 2) . The strongest ten weighting OTUs for each surface type are 
highlighted if they were also significant indicator OTUs (all lvalues <0.05). (c) All samples were compared to potential source environments using 
principal components analysis (PCA), and the first principal component (37.8% of variance explained) was used as a surrogate for community simi- 
larity to either phyllosphere or human skin bacterial communities. Boxplots delineate (from bottom) minimum value, Ql, median (Q2), Q3, max- 
imum value; notches approximate 95% confidence around median value, and outliers fall outside of the quartile range. Letters above each box 
indicate significant groupings after Tukey's hones significant difference (HSD) test (adjusted lvalue <0.05). 
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degrees of influence, or weight, on DB-RDA clustering 
patterns, and we selectively tested the strength of associ- 
ation between the strongest weighted individual bacterial 
OTUs and relative surface types (that is, indicator value, 
as defined in [20]) with the indval function in labdsv. 
Mantel tests, using the mantel function in vegan, were used 
to test for a community dissimilarity correlation with spatial 
distance within the classroom for each of the four surface 
types, as well as all sites regardless of surface type. To assess 
community similarity to various potential source habitats 
(soil, aquatic, phyllosphere, human skin and human 
gut), we compared all samples at the class level to 
publicly-available datasets from representative potential 
source environments [4,21-24] using principal compo- 
nents analysis. Positions of each point in relation to source 
environments along the first principal component are 
used to compare their similarity to either human skin or 
phyllosphere bacterial communities. 

Results and discussion 

After rarefaction to 4,000 sequences per sample, 58 sur- 
face samples were represented by a total of 3,745 bacterial 



OTUs. Approximately half (51.4%) of OTUs were found 
only once or twice, and these were included in analysis. 
The most commonly detected OTU {Sphingomonas sp.) 
accounted for 1.9% of all sequences. The four surface types 
harbored significantly different communities {P - 0.001; 
F = 2.34; from permutational multivariate analysis of 
variation on Canberra distances; Figure 2). Within each 
surface type, bacterial community composition was not 
significantly predicted by spatial distance among samples 
(all four P values >0.1; from Mantel tests). Even when all 
surfaces were considered together, spatial distance was not 
a significant predictor of bacterial community composition 
(P = 0.4). These results suggest that site-specific factors (for 
example, habitat selection or dispersal from specific body 
sites) are more important than dispersal among sites for 
bacterial community assembly in the built environment. 

We found that surface communities were significantly 
associated with taxa specific to distinct microbial sources 
(Table 1). Bacterial OTUs labeled in Figure 2b are the 
most strongly weighted OTUs for each surface type that 
are also significant {P < 0.05) indicator OTUs. Three of the 
strongest indicators for chair seats were Lactobacillus 



Table 1 Closest known isolates related to indicator operational taxonomic units 



Greengenes genus 


P value 


Surface type 


Closest 16S NCBI isolate 
and accession 


Isolate source 
environment 


Sequence similarity 
to isolate (%) 


Lactobacillus 


0.00 r 


Chairs 


Lactobacillus John son ii NR_075 064.1 


Human gut 


99 


Corynebacterium 


0.001* 


Chairs 


Corynebacterium resistens NR_040999.1 


Human infection 


99 


Corynebacterium 


0.00 r 


Chairs 


Corynebacterium confusum NR_026449.1 


Human clinical specimens 


99 


Staphylococcus 


0.01 r 


Chairs 


Staphylococcus epidermidis NR_0 74995.1 


Human skin 


99 


Corynebacterium 


0.001* 


Chairs 


Corynebacterium riegelii NR_026434.1 


Human urinary tract 


99 


Staphylococcus 


0.019* 


Chairs 


Staphylococcus saprophyticus NR_074999.1 


Human urinary tract 


99 


Lactobacillus 


0.00 r 


Chairs 


Lactobacillus crispatus NR_074986.1 


Human vagina 


99 


Lactobacillus 


0.003* 


Chairs 


Lactobacillus acidophilus NR_075049.1 


Human gut 


99 


Streptococcus 


0.001* 


Desks 


Streptococcus oralis NR_1 02809.1 


Human oral 


99 


Streptococcus 


0.001* 


Desks 


Streptococcus salivarius NR_1 0281 6.1 


Human oral 


99 


Brevundimonas 


0.002* 


Desks 


Brevundimonas variabilis NR_037 106.1 


Pond water 


99 


Streptococcus 


0.001* 


Desks 


Streptococcus intermedius NR_1 02797.1 


Human purulent infection 


99 


CandidatusPhytoplasma 


0.001* 


Desks 


None** 






Alicyclobacillus 


0.001* 


Walls 


Tumebacillus permanentifrigoris NR_043849.1 


Soil 


99 


Chroococcidiopsis 


0.028* 


Walls 


Halospirulina tapeticola NR_026510.1 


Saline aquatic 


96 


Alicyclobacillus 


0.001* 


Walls 


Tumebacillus permanentifrigoris NR_043849.1 


Soil 


98 


Rhodopseudomonas 


0.001* 


Walls 


Methylobacterium adhaesivum NR_042409 


Drinking water 


98 


Salmonella 


0.001* 


Floors 


Pantoea ananatis NR_1 03927.1 


Phyllosphere 


99 


Roseomonas 


0.001* 


Floors 


Roseomonas gilardii NR_029061.1 


Human blood 


99 


Roseomonas 


0.001* 


Floors 


Roseomonas frigidaquae NR_044455.1 


Water-cooling system 


99 


Salmonella 


0.001* 


Floors 


Pantoea ananatis NR_1 03927.1 


Phyllosphere 


99 



All extant operational taxonomic units labeled in Figure 2 (and thus influential in distance-based redundancy analysis, as well as significant indicator taxa for their 
respective surface type) were related to their closest known bacterial isolate using 16S rRNA sequences in the NCBI Bacteria & Archaea Isolate Database. Source 
environments are from each isolate's respective published source environment. ^Unadjusted P value < 0.05. ^^Closest known isolate 89% similar. NCBI: National 
Center for Biotechnology Information. 
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OTUs closely related to taxa commonly found in the hu- 
man gut and vagina, in addition to other human-associated 
Corynebacterium and Staphylococcus species. Desk surfaces 
were significantly associated with several indicator taxa 
from human body habitats; two of the strongest desk sur- 
face indicator OTUs were Streptococcus species commonly 
found in human skin and oral samples, as well as another 
Streptococcus originally isolated from purulent infections. 
Floors in the classroom also harbored skin-associated 
OTUs, but were more strongly indicated by a cyanobacter- 
ial OTU, presumably from a non-human environmental 
source. Walls are likely the surfaces in this study that have 
the least contact with humans; this was consistent with our 
observation that these surfaces were associated with indi- 
cator taxa related to Sphingomonas and Alicyclobacillus 
species that are commonly abundant in airborne bacterial 
assemblages [10,25]. Both walls and floors held larger rela- 
tive proportions of Cyanobacteria than seats and desktops, 
likely reflecting soil and bioaerosols as microbial sources. 
Few other major phylogenetic differences are evident at a 
broad taxonomic level (Figure 3). 



The community variation we observed among surface 
types could be the result of either surface- specific envir- 
onmental filtering or dispersal. Evidence to date suggests 
that microbial communities on indoor surfaces, with a few 
exceptions, are likely primarily dispersal-driven given the 
well-documented inputs from dust, settled airborne par- 
ticles, outdoor air and surface contact [7-9,25-27]. The 
relative importance of in situ community dynamics and 
habitat filtering remain to be fully understood in the 
built environment. 

Dispersal between human skin and contacted surfaces 
is well documented, primarily in a medical context using 
known pathogenic strains [1,28,29]. However, the vast 
majority of bacteria on and in the human body are not 
pathogens, but rather appear to be commensal components 
of our own microbiome [30]. Recent evidence has sug- 
gested that the unique microbial assemblages detected on 
different human body parts can be transferred to indoor 
surfaces following contact [26,31], and that these bacterial 
traces of human contact are evident in places such as 
restrooms [7]. Our results suggest that this transmission is 
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Figure 3 Taxonomic composition of all 58 samples used in this study. Samples are grouped by surface type. All taxonomic groups representing 
<1% (Phylum and Class) and <4% (Order) of sequences were grouped into 'other'. 
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possible even when that contact is indirect, and that such 
contact has a greater impact on microbial community 
structure indoors than dispersal among surfaces. 

Conclusion 

We are constantly surrounded by diverse microbial com- 
munities indoors, and we are just beginning to understand 
how our interactions in the built environment shape those 
communities and our own human microbiome. Our results 
indicate that human-associated microbial communities can 
be transferred to indoor surfaces following contact, and that 
such transmission is possible even when contact is indirect, 
but that proximity to other surfaces in the classroom does 
not influence community composition. 

Availability of supporting data 

Sequence files and metadata for all samples used in this 
study have been deposited in Figshare (http://dx.doi.org/ 
10.6084/m9.figshare.687155). A flail record of all statistical 
analysis is included as Additional file 5, and was created 
using the knitr package in R [19]. Original R scripts are 
available in GitHub (https://github.com/jfmeadow/Mea- 
dow_etal_Surfaces). Metadata, the unrarefied OTU table, 
and corresponding taxonomic classifications have all been 
included as Additional files 2, 3 and 4, respectively. 
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